Structural signatures of the unjamming transition at zero temperature 
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We study the pair correlation function g(r) for zero-temperature, disordered, soft-sphere packings 
just above the onset of jamming. We find distinct signatures of the transition in both the first and 
split second peaks of this function. As the transition is approached from the jammed side (at higher 
packing fraction) the first peak diverges and narrows on the small-r side to a delta-function. On 
the high-r side of this peak, g(f) decays as a power- law. In the split second peak, the two subpeaks 
are both singular at the transition, with power-law behavior on their low-r sides and step-function 
drop-offs on their high-r sides. These singularities at the transition are reminiscent of empirical 
criteria that have previously been used to distinguish glassy structures from liquid ones. 

PACS numbers: 61.43.Fs, 81.05.Rm, 64.70.Pf, 



It is only natural to suspect that the dramatic dynam- 
ical arrest that occurs as a liquid is cooled into a glass 
must be accompanied by a signature in the underlying 
atomic arrangements. However, the atomic configura- 
tions in the liquid and glass are strikingly similar to one 
another. Over the years, the challenge to identify a sub- 
tle structural difference between the two states has led to 
the proposal of several empirical criteria 0, . In this 
paper, we revisit an old idea, due to Bernal jj], of using 
static sphere packings to gain insight into the structure 
of amorphous systems [a, B S B S Li_LJ| ■ We find that, 
with decreasing density, the structure of such packings 
changes distinctly as they unjam. 

It is particularly revealing to study static packings of 
spheres interacting via a finite-ranged, purely-repulsive 
potential. There is a fundamental change in the me- 
chanical properties in such systems, reminiscent of the 
glass transition as the packing fraction <j> is varied 
across a well-defined unjamming/jamming transition at 
4> c , which was found to coincide with the value of ran- 
dom close packing 0.64 [lj, Above <f> c the sys- 
tem has nonzero static shear and bulk moduli, while be- 
low <f> c it costs no energy to shear or compress the sys- 
tem by an infinitesimal amount. Moreover, the unjam- 
ming transition in many ways resembles a critical point 

with many quantities, 
including a diverging length scale, behaving as power- 
laws as the transition is approached from higher density. 
Here we show that the unjamming transition exhibits 
clearly identifiable structural signatures associated with 
diverging quantities, even though both the jammed and 
the unjammed states are disordered. These structural 
characteristics are echoed in some of the empirical crite- 
ria P, 0, that have previously been proposed for the 
glass transition. 

The simulations reported here are for monodisperse, 



soft spheres of diameter a that interact through the po- 
tential V[r) = (Vo/a)(l - r) a , for r < 1, and V(r) = 
when r > 1. Here, r is the center-to-center separation 
between two particles, measured in units of a. We have 
studied both the harmonic, a = 2, and the Hertzian, 
a = 2.5, cases. Particles are defined to be in contact if 
they overlap. Our three-dimensional systems consist of 
1024 < N < 10000 spheres in periodic, cubic, simulation 
cells. To enable a systematic study of the approach to the 
unjamming transition, we employ conjugate-gradient en- 
ergy minimization |21| to obtain T = configurations at 
various packing fractions <j>. We average over ensembles 
of configurations at the same distance from the transi- 
tion point, i.e., at the same values of A</> = (j)— <f> c , which 
is equivalent to averaging over systems with the same 
pressure [l2llT3| . 

Structural signatures of jamming are more evident 01 
in the pair correlation function, <?(r), than in the struc- 
ture factor, S(k), even though the two functions are sim- 
ply related by a Fourier transform. We therefore focus on 
the structural characteristics that signal the approach of 
the zero-temperature transition from the jammed side, 
through a detailed analysis of g(r) computed with ex- 
traordinary resolution. In particular, we study two fea- 
tures of g(r) that undergo dramatic changes at the un- 
jamming transition: the first peak and the split second 
peak (z. e., the two subpeaks that merge with increas- 
ing temperature to become the second peak in a typical 
dense liquid). A plot of g(r) is shown in Fig. ^ above the 
jamming transition. 



I. FIRST PEAK OF g(r) 

The dominant feature in g(r) is the first tall, sharp 
peak at r pca k (Fig. [IJ. Precisely at the jamming thresh- 
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FIG. 1: The pair correlation function g(r) vs. r at two ex- 
treme values of Acf> = <f> — <f> c — 10~ 6 (dotted line) and 10~ 
(solid line). The maximum value of the first peak height is 
higher and its width narrower for the lower value of Ac/>. For 
Acj> — 10 -6 , the first peak maximum is approximately 10 6 , 
far beyond the scale of the graph. 



old, Acf> = </) — <j) c = 0, this peak is a <5-function at 
^pcak = 1; g{r) is precisely zero for r < 1 and has a 
power-law tail extending to r > 1. The weight under 
the <5-function is the coordination number at contact, 
^contact- (As we discuss below, Z contac t is a few percent 
less than the isostatic coordination number, Z„ = 2d = 6 
for our d = 3-dimensional systems EU E2-) For 
A<fi > 0, there is some overlap between particles so that 
the delta- function peak broadens and shifts to r pca k < 1 . 
The broadening produces a tail extending to r < r-p Ca k 
that disappears in the limit where Acfi vanishes. 

By analyzing the height of the first peak and its left- 
hand width, we showed |l3j | that the peak approaches a 
delta-function as Acf> decreases towards zero. We have 
since obtained more systems over a wider range of A<f>. 
Figure [21 shows the dependence on Acfi of the first-peak 
height (/(fpeak) and its left-hand width wl for both har- 
monic and Hertzian potentials. Independent of the inter- 
action potential, we find: 



^(fpcak) ~ •' 



(1) 

(2) 



These scalings are consistent with the area of the peak 
approaching a constant in the limit Acj) — > 0. 

We turn now to the shape of the first peak in g(r) 
at r < 1 for a system at A<fi = 1 x 10~ 6 , just above 
the jamming transition. As shown in Fig. |2Ia), for the 
region r < 1, g(r) is almost strictly exponential with only 
a small curvature near its peak. It can be fit with the 
functional form: 



g(r < a) = g exp - 
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FIG. 2: The height, <?(r pca k), (top panel) and the left-hand 
width, wl, (bottom panel) of the nearest-neighbor peak of 
g(r), over several orders of magnitude of Acj), for monodisperse 
spheres with purely repulsive (a) harmonic spring and (b) 
Hertzian interactions. 



As we compress the system above <f> c we see that the 
exponential behavior of the tail at r < 1 gradually be- 
comes more Gaussian as the system is compressed above 
the transition. This is shown in Fig. |3^-d. We can still 
use Eq. [21 to fit the shape, however, with different coeffi- 
cients a.\ and 0.2- In Fig. 21 we show the evolution of ot\ 
and «2 with A</>. This evolution occurs more rapidly for 
Hertzian (not shown) than for harmonic spheres. 

The results for g(r < 1) have implications for the dis- 
tribution of inter-particle normal forces, P(F). This 
is shown in Fig. At high compressions, P(F) is 
well-described by a Gaussian, but the tail straightens 
out towards an exponential as A<p is lowered towards 
zero. These results are consistent with previous results 
of Makse, et al. 0, |2^|, who studied sphere packings 
at fixed pressure. (As we noted above, constant pres- 
sure corresponds to constant A</>.) The Gaussian shape 
at high Ac/) is consistent with expectations for equilib- 
rium systems interacting with a harmonic potential |24j . 
However, these systems are at zero temperature, and it is 
unclear whether they can be described by a nonzero effec- 
tive temperature. The exponential behavior at small A</> 
agrees with experimental and simulation data on static 
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FIG. 3: Pair distribution function g(r) up to contact, < r < 
1, on a linear-log scale, for A4> = (a) 10~ 6 , (b) 1CT 4 , (c) 1CT 2 , 
and (d) 10" 1 . 
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FIG. 4: Evolution of the parameters ai and «2 (top) and g 
(bottom) in Eq. with A(j>, for harmonic repulsions. The 
parameters 01,2 are clearly related to the pressure, while g a 
is associated with g(r pea k). 



granular packings of hard partic les, which necessarily ex- 
ist at packing fractions near (j> c 0, [2f|, [2(J . 

There is interesting behavior above the asymmetric 
first peak in g(r) as well as below it. Fig. HJa) shows 
that g(r > 1) versus r — 1 varies as a power-law for a 
system just above the transition at A<fi = 10~ 8 : 



l{r > 1) oc [r — 1] 



(4) 



with rj = 0.48 ± 0.03. This result was first reported for 
gravity-sedimented, granular packings [2]j, but over a 
much smaller range in g(r) than presented here. We note 
that there is a very slight knee that occurs near r — 1 = 
3 x 10~ 2 . The asymptotic power law behavior near r = 1 
should be determined only from the region below this 
knee. As we will show below, this knee becomes more 
pronounced as A<p increases. 

The number of neighbors, Z{€), that are separated by 
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FIG. 5: The distribution of normal contact forces, P(F), for 
the purely repulsive, harmonic potential, at different compres- 
sions A<j>. 



a distance of at most i [28| is given by the integral 



Z(£) = 240 / g{r')r' dr' . 



(5) 



This is shown in Fig.Efb). Therefore, Fig.E^a) and Eq.0] 

imply, that for a system at the transition, Z should in- 
crease with I as 



Z — Z c , 



^l-r f =0.52±0.03 



(6) 



where r\ is defined in Eq. 0] and Z CODtac t — 5.88 is the 
average number of neighbors per particle at contact at 
the transition. This scaling is consistent with the one 
reported by O'Hern et al. who looked at how the 

excess coordination number increased as a system was 
incrementally compressed above (f> c . We note that we 
have found a similar exponent of 0.50 ± 0.03 using the 
Hertzian interaction potential, a — 2.5, as we found for 
the harmonic potential, a = 2. 

Donev et al. |2£| also calculated the number of neigh- 
bors Z(€) near the transition, but for a system of hard 
spheres. When they plotted \o%{Z{€) — Z c ) versus log£ 
they found that the slope was closer to 0.6 than 0.5. They 
argued that this implies that g(r > 1) oc (1 — r)~ QA for 
hard spheres, which is different from our result for soft 
repulsive spheres (Eq.|BJ. Based on this difference, they 
concluded that the power-law behavior seen for r > r pca k 
is not universal and depends on the interparticle poten- 
tial as well as perhaps the procedure leading to the T = 
configuration. 

However, that conclusion does not follow from their 
analysis. The confusion occurs because their results were 
for a system in which all the rattlers (that is, particles 
without any neighbors) had been removed. It is not sur- 
prising that this changes the nature of the pair corre- 
lation function from a system, such as ours, where all 
particles are considered. That is, we analyze the system 
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FIG. 6: Behavior just above the first-neighbor peak for a 
system at A(f> = 10 -8 . (a) g(r > I) versus r— 1. A power law 
with exponent of —0.5 is indicated, (b) Z(£) vs. i = r — 1, 
computed by numerically integrating g(r) > 1. A power law 
with exponent of 0.5 is indicated. 



directly as produced by energy minimization without per- 
forming the additional step of removing particles that do 
not happen to be part of the backbone of the structure. 
Whether one chooses to remove rattlers or not depends 
on the physics one wants to study. Here we remark that 
it is more natural to include rattlers in g(r) if one is going 
to compare to any experiment. We also note that rattlers 
can join the backbone when the system is perturbed, and 
can therefore influence the systems response and stabil- 
ity, as argued in Ref. [30j . Finally, we claim that the 
conclusion of Donev, et al. that the power-law exponent 
depends on the potential of interaction or the algorithm 
for creating the states (aside from removing one class of 
particles by hand) is unfounded. Rather, the difference 
arises from whether one studies systems with rattlers in- 
cluded or removed. When Donev, et al. include rattlers 
in their analysis, they also find an exponent closer to 
0.5 |3lj . consistent with what we have found here and in 
previous work [l^, EH H3 ■ 

There is a difference between Z c and Contact- Z c is 
the number of contacts only for a system that has all 



rattlers removed from the system. If the rattlers are not 
removed from the system then Z contac t < Z c . Because 
g(r), for example measured in an experiment, does not 
distinguish between rattlers and particles belonging to 
the connected backbone of the structure, one must use 
^contact ~ 5.88 in Eqs. El Indeed, calculating Z{t) either 
directly as in Ref. [2!j, or integrating g(r) via Eq. \5\ are 
identical only when using Z contac t = 5.88, and not Z c = 6 
in Eq. El If N denotes the total number of particles, 
^overlaps the number of overlapping pairs, and Nf the 
total number of rattlers, then Contact = 2./V ovcr i aps /7V 
while Z c = 2N ovclUps /(N - Nf). Although only 2% of 
the particles are rattlers at the transition, their exclusion 
can produce a change of the power law exponent. 

We find that the observed power-law for g(r) depends 
on A(p (see Fig. 0a)). As Atfi increases, the knee in 
g(r) near r — 1 ~ 3 x 10~ 2 becomes more prominent. 
In the region below this knee (i. e., at smaller r), the 
slope of log[<7(r)] versus log[r — I] decreases. We show 
this slope, 7], as a function of A<fi in Fig. 0b). As A(j> 
approaches zero, the value of r\ increases and approaches 
i] = 0.5. A similar trend was noted in the X-ray tomog- 
raphy experiments of Aste et al. j32, who measured g(r) 
inside large, 3D granular packings. Note that Donev, et 
al. p9l | study configurations below the unjamming tran- 
sition, and it is not clear if the apparent value of rj also 
changes as the density is decreased below C . Fig. [7{c) 
shows, for different values of A(p, how g(r) behaves in 
the region from r = r poa k out to the first minimum at 
r ~ 1.4. There are several notable features in g(r) in the 
vicinity of contact that are apparent when g(r) versus 
(r — ^poak) is plotted on log — log axes. There is a drop 
in g(r) that occurs at r — 1 for each value of A(j> (note 
that this corresponds to different values of r — 7' pea k for 
different Aip). The separation r = I distinguishes parti- 
cles that are overlapping from those which are just out 
of contact. The magnitude of the jump decreases, and 
the extent of the power-law region described by Eq. 0] 
also decreases, as the system is progressively compressed 
above the unjamming transition. The region beyond con- 
tact is relatively unaffected by compression. (Although 
the power law exponent changes slowly with A<f> as high- 
lighted in Fig. EJb), by far the largest change occurs in 
crossing from r < 1 to r > 1.) This indicates that as the 
system is compressed, particles are depleted from the re- 
gion beyond contact, r > 1, and are absorbed into the 
contact region, r < 1. 

We have determined that this drop in g(r) is not an 
artifact of the zero-temperature system. We have used 
several different protocols for creating particle configu- 
rations, as shown in Fig. |S| and described in its caption. 
Evidently, the drop in g(r) persists to finite temperature. 

Why should there be an abrupt drop in g(r) at r = I? 
One possibility is that, to a first approximation, compres- 
sion only changes the nearest-neighbor spacing of parti- 
cles that are already in contact, that is, within a distance 
r < 1 of each other. Particles separated by a distance 
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FIG. 7: Pair distribution function g(r), just above the first 
peak for A(j> = 10" 8 (x), I0~ 6 (circle), I0" 4 (squares), 
10~ 3 (diamonds), 10~ 2 (+). (a) g(r) vs. r — 1, showing the 
power-law behavior in (r — I), (b) The exponent r\ character- 
izing the power-law behavior of (a), as a function of A<j>. (c) 
g(r) from the nearest-neighbor, first peak at r pea k out to the 
first minimum at r « 1.4 as a function of r — r pca k. 
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FIG. 8: Pair distribution function g(r), from the nearest- 
neighbor, first peak at r- pca k out to the first minimum at 
r fa 1.4, at A(j> = 10 -4 . Different symbols represent different 
configuration protocols with and without temperature: conju- 
gate gradient minimization to T — (open circles), quenched 
molecular dynamics at T — (open squares), molecular dy- 
namics at a very low temperature (solid diamonds) , molecular 
dynamics at a higher temperature (+). The jump in g(r) per- 
sists to small but nonzero temperature, and is smoothed out 
at high enough temperatures. 



lar to the line connecting them. This is consistent with 
our data. 

As the system is compressed and particles are incorpo- 
rated into the contact region, the number of overlapping 
pairs increases. As mentioned above, the average number 
of overlapping neighbors per particle, Z , increases with 
compression as Z — Z contac t oc A0 5 . As this occurs, the 
distribution of Z values, P(Z), also shifts. This is shown 
in Fig. |5J Close to the transition, most particles have 6 
overlapping neighbors. As the system is compressed to 
A<p = 0.1, the maximum shifts to Z = 9 but there is still 
no observable weight at Z = 12. 
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r > 1 may not be moved appreciably closer to one an- 
other by compression. Instead, we suggest that upon 
compression, the movement of particles that are not yet 
overlapping is predominantly in a direction perpendicu- 



FIG. 9: The distribution of the number of overlapping neigh- 
bors per particle, Z, at several different compressions. Key 
same as Fig. El 
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II. SPLIT SECOND PEAK OF g(r) 

The approach of the unjamming point is also apparent 
in features associated with the second peak in g(r). Fig- 
ure El shows this region of g{r) for several values of A(f>. 
For the system closest to jamming, at A0 = 1 x 10~ 6 , 
fFig. HUf a')') there is a pronounced splitting of the second 
peak in g(r) into two sub-peaks, located at r = y/3 and 
r = 2. Such a splitting of the second peak has long been 
known Indeed, the emergence of a split second peak 
was one of the early criteria used to signal the onset of 
the glass phase in supercooled liquids 2] . 
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FIG. 10: g(r) region around the split second peak for A(j> = 
(a) 10~ 6 , (b) 10" 3 , and (c) 10 _1 . 

On the small-r side of both of these sub-peaks we find 
that g(r) increases rapidly. In Fig. El we have attempted 
to characterize these features for the system closest to the 
transition. In that figure we plot g(r) versus log(v3 — r) 
and log(2 — r). We have not been able to tell whether g(r) 
itself diverges or whether the divergence appears only 
in its slope, dg(r)/dr. We would need to average over 
many more systems to tell these two possibilities apart 
unambiguously. Our best fits to the two possible cases 
are: 

g(r) = ai (V3-r)- Q2 (7) 
g(r) = g(V3)-h(V3-r) h > (8) 
g(r) = - Cl log(V3-r) (9) 

g(r) = d x (2-r)-* (10) 
g(r) = g(2)- ei (2-r) e2 (11) 
g(r) = -/!log(2-r) (12) 

where the fit parameters are provided in Table [I] 



and 



TABLE I: Fit parameters for the two sub-peaks, situated at 
r = v3 and r = 2, respectively, that make up the split second 
peak in g(r). For = 10" 6 . 
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ai = 0.9 


di = 1.1 
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FIG. 11: Functional form of the sub-peaks of the split-second 
peak in g(r) for Atf) = 10 -6 . Left-hand sub-peak at r = v3- 
Right-hand sub-peak at r — 2. The solid line corresponds to 
a power-law fit, the dashed line to a logarithmic fit, and the 
grey line to a shifted power- law fit. 



On the high-r side of each sub-peak, it is clear from 
Fig. HOfa) that there is both a step-function drop-off 
[29l l33| and an additional smooth decrease of g(r) with 
increasing r. Although, from these data, we cannot un- 
ambiguously determine a functional form for this last 
smooth decrease, the data suggest that it might be fit 
with 



g(r) = a - (r - r s ) b 



(13) 



where < b < 1 and r s — V3 and 2 for the first and 
second sub-peaks, respectively. As the system is com- 
pressed above the transition, the structure around both 
sub-peaks becomes rounded. In particular, it is only in 
the limit of A0 — > that the drop-off becomes a sharp 
step- function. 

Here, we suggest a simple interpretation of the ori- 
gin of the step- function drop-off in g(r) on the high-r 
side of both sub-peaks. At the jamming transition, each 
particle must be held in place by particles that are just 
in contact with no overlap. At that point there is an 
average of precisely Z c = 6 neighbors per particle in 
the force network. It seems plausible that the second- 
nearest-neighbor peaks originate from pairs of particles 
that have at least one neighbor in common, while pairs 
that do not have a common neighbor contribute only to 
a slowly-varying background. 

The separation r = 2 corresponds to two particles 
on diametrically opposite sides of a common neighbor 
(Fig. El right figure) . This is the largest distance that 
can separate two particles that have one common neigh- 
bor. The separation r = V3 corresponds to the largest 
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possible separation between two particles, i and j, that 
have two common neighbors (Fig. El left figure). In 
this case the four particles lie in the same plane; the two 
common neighbors touch and the particles i and j are 
on opposite sides of the crease between them. When r 
is only slightly greater than 2, there can no longer be 
any contribution from two particles that share a com- 
mon neighbor. Similarly, when r is only slightly greater 
than a/3, there can be no contribution from pairs of par- 
ticles that share two common neighbors. This leads to 
step- function drop-offs at r = a/3 and r — 2. 



ccc 



FIG. 12: Left: The particles at the left and at the right share 
two common neighbors and are separated by a distance r — 
a/3. Right: The particles at the left and the right share one 
common neighbor and are separated by r = 2. 

Support for the hypothesis that the drop-offs are due 
to pairs of particles that share common neighbors can 
be obtained from the angular correlations between two 
neighbors of a common particle. We show this angular 
correlation, P(9), in Fig. El where 9 is defined in the 
sketch in the inset. At 6 = tt (which corresponds to 
particles separated by r = 2) Fig. I13f d) shows that the 
distribution goes to zero as, 



P(9)~(n 



\x=0.75 



(14) 



If the sub-peak at r — 2 arises from pairs that share one 
common neighbor, then the form of g(r) just below r = 2 
should be related to the form of P(9) just below 6 = tt 
by 



g (r) = P(9)d9/dr 
where the Jacobian factor diverges as 



dd^j/dr 



Thus, Eg. H*5l implies 



1 



cos -|i a/4 - r 2 



g(r) ~(2-r) ; 



(15) 



(16) 



(17) 



This is in reasonable agreement with the fit in Eq. 1111 
where we found g(r) ~ (2 — r )~ d 2=-o.08^ j m piy m g that 
the second sub-peak does indeed arise from two particles 
that share a common neighbor and that this leads to the 
observed step-function drop-off just above r = 2. How- 
ever, we should note that this argument does not neces- 
sarily imply that g(r) diverges at r = 2, as indicated by 



Eq. 1171 since we cannot determine the behavior of P(9) 
asymptotically close to n with sufficient accuracy. 

By a similar argument, if the first sub-peak arises from 
pairs of particles that share two common neighbors, then 
the peak at r = \/3 corresponds to 9 — 2tt/3. For 9 < 
2ir/3, we find 



P{9) ~ (2tt/3 - 6>) 



-0.17 



(18) 



This behavior implies g(r) ~ (Vo — r)~ 017 ; this result 
is not too different from the fit in Eq. 7, which suggests 
g(r) ~ (a/3 — r)~ ' 12 . This consistency check suggests 
that the first sub-peak does indeed arise from pairs of 
particles that share two common neighbors. 
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FIG. 13: (a) The three-particle angular correlation, P(6i,j), 
versus Gij for A(f> = 10 -6 . Bottom four panels are fits to the 
regions around; (b) 9 = 60°, (c) 6 = 120°-left, (d) 9 = 120°- 
right, (e) 9 = 180°. The power law exponents are ~ 0.5, 0.17, 
0.17, and 0.75, respectively 

Although pairs with two common neighbors cannot 
contribute to g(r) for r > a/3, pairs with only one com- 
mon neighbor can still contribute. Such pairs may ac- 
count for the smooth decrease in g(r) (see Eq. I13|) just 
above r = a/3. 
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III. DISCUSSION 

We have shown that the unjamming transition is ac- 
companied by several features in the pair correlation 
function: 

• A delta- function at r — 1 with area ^contact < 6. 

• A power-law at r > 1 of the form g(r) ~ (r— 1) — °" 5 . 

• A sub-peak at r = \/3 that either diverges or has 
diverging slope as r — + V3 and that has a step- 
function drop-off just above r = y/3. 

• A sub-peak at r — 2 that either diverges or has 
diverging slope as r — * 2~ and that has a step- 
function drop-off just above r = 2. 

These features appear for both harmonic and Hertzian 
repulsions, and therefore seem to be purely geometrical 
features of the jamming/unjamming transition at zero 
temperature. 

Here, we suggest that some of these structural features 
are reflected in a less extreme form in the two primary 
empirical criteria that have been used extensively in the 
literature to identify the glass transition. The first is the 
ratio, TZ, of the first minimum to the first maximum in 
g(r). As the temperature T of a glass- forming liquid de- 
creases, TZ decreases. Wendt and Abraham [2] proposed 



that TZ = 0.14 corresponds to a reasonable estimate of 
the glass transition temperature, T g . We note that at 
the unjamming transition we have studied here, TZ = 0. 
It is possible that the decrease in TZ observed as T is low- 
ered towards T g is a remnant of the vanishing of TZ that 
occurs at the unjamming transition. The latter prop- 
erty may be the underlying reason for the success of this 
empirical criterion. 

A second popular empirical criterion concerns the sec- 
ond peak of g(r). As T is lowered towards T g , the second 
peak splits into two sub-peaks. The flattening of the sec- 
ond peak that signals its splitting has been used to iden- 
tify T g . We suggest that the splitting may reflect the 
singular sub-peaks that occur at jamming/unjamming 
transition. The singular nature of the splitting at this 
transition may provide the fundamental underpinnings 
of this criterion. 
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